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^ Abstract 

We use electric dipole spin resonance to measure dynamic nuclear polarization in InAs nanowire 
quantum dots. The resonance shifts in frequency when the system transitions between metastable 
high and low current states, indicating the presence of nuclear polarization. We propose that the low 
and the high current states correspond to different total Zeeman energy gradients between the two 
quantum dots. In the low current state, dynamic nuclear polarization efficiently compensates the 
Zeeman gradient due to the (/-factor mismatch, resulting in a suppressed total Zeeman gradient. 
We present a theoretical model of electron-nuclear feedback that demonstrates a fixed point in 

o 

nuclear polarization for nearly equal Zeeman splittings in the two dots and predicts a narrowed 
hyperfine gradient distribution. 
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Hyperfine interaction couples electron spin to nuclear spins enclosed by the electron's 
wave function. In the context of spin qubits in III-V semiconductors, the most prominent 
effect of this interaction is that fluctuating nuclear spins cause electron spin dephasing [IH3] • 
Interestingly, ideas for suppressing nuclear spin fluctuations also rely on the same hyperfine 
interaction, since electron spin transport can lead to dynamical nuclear polarization (DNP) 
in quantum dots [3]. An experimental manifestation of DNP is a hysteretic current in the 
spin blockade regime |5H8] . Sometimes, including in the present work, the hysteresis could 
be extended to high magnetic fields, suggestive of a large degree of nuclear polarization 

pun]. 

In this report, we study hysteretic spin blockade in InAs nanowire quantum dots using 
electric dipole spin resonance (EDSR) spectroscopy [T1HT3] . Surprisingly, the degree of 
polarization deduced from EDSR does not exceed a few milliTesla, much smaller than the 
hysteresis range. We explain this apparent contradiction by nuclear fields compensating the 
natural Zeeman energy difference between the two quantum dots caused by the mismatch 
of their ^-factors [2|. In this case one of the double dot states (T ) is blocked leading to 
a reduced current. We support this idea by analytical and numerical calculations of spin 
blockade transport in the presence of hyperfine and spin-orbit interactions. 

From our model we deduce a narrowing of the hyperfine gradient distribution to a few 
percent of the unpumped distribution width. This finding is especially relevant for two- 
electron singlet-triplet qubits, where the hyperfine gradient is the source of dephasing [15J. 
We predict an order of magnitude enhancement in the coherence time due to gradient sup- 
pression induced by spin blockade transport. This is an alternative route to T 2 * enhancement 
compared with nuclear spin pumping by pulsing the double dot through an S-T transition 

pMEE]. 

We present data from two devices that were studied in two previous publications [3j [19] . 
InAs nanowires, 40-80 nm in diameter, are deposited on top of five narrow bottom gates 
which are then used to define few electron double quantum dots. The nanowires are con- 
tacted by Ti/Al leads to measure electron transport through the system. All measurements 
are performed at T = 250-300 mK in 3 He refrigerators. 

We tune the double dots to the so-called strong coupled spin blockade regime, in which 
the current exhibits a dip rather than a peak at zero magnetic field due to an interplay 
between tunnel coupling, hyperfine and spin-orbit interactions [TUl IT9"] . Fig. 1(a) shows the 
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FIG. 1: (a) Charge stability diagram near the (3,1)— >(2,2) transition (Vd c = 7 mV, B ext = 0.79 
T). Low and high current states are labeled, the dashed arrow indicates the detuning axis e. (b) 
Dependence of the leakage current on magnetic field and detuning. The detuning is swept along 
the line in panel (a), and the field is stepped after each sweep. For (a) and (b) the sweep direction 
is from left to right, (c) Magnetic field retrace at fixed detuning e ~ 3.5 meV, the sweep directions 
are indicated with arrows. The current becomes zero above 4 T since the double dot then shifts 
into a Coulomb blocked state due to the Zeeman shift of T + (l,l). The same device was previously 
studied in [T9] . 

double dot current in the vicinity of a spin-blocked charge degeneracy point, which appears 
as a double-triangle shape when the left and right gates are swept at a finite d.c. voltage 
bias Vdc across the double dot. Inside the triangles, one can see sudden transitions between 
a low and a high current state. The difference in current between the two states exceeds, 
for some settings, an order of magnitude (1-2 pA vs. 10-20 pA). This switching is only 
visible within the boundaries of the triangles. The boundaries themselves remain fixed in 
gate voltage. This indicates that the jumps in current do not originate from charge switches 
in the vicinity of the double dot. 

The appearance of the two current states is strongly influenced by an applied magnetic 
field -Bext [Fig. 1(b)]. At zero field we always observe a stable low current because we are 
in a zero-field dip characteristic of all strongly coupled quantum dots [5j [191 120]. In the 
absence of hysteresis this dip is hundreds of mT wide with current increasing smoothly with 
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field. However, in the regime where we observed the switching, the current drops when the 
field increases beyond 10-20 mT, and the system enters a metastable low current state [1 0J . 
At higher fields, both low and high current states are observed, as well as sharp transitions 
between the two. 

In Fig. 1(c) we present an example of a hysteretic current trace. We see that the low 
current state can be "dragged" up to very high fields (over 4 T). The reversed sweep shows 
a distinct high current state down to ~ 0.7 T, where the current switches. If the magnetic 
field is fixed to a value inside the hysteretic regime, the double dot may suddenly switch after 
minutes or seconds, or remain in either of the states for as long as hours (see supplementary 
material for time-dependent measurements). We note that the higher noise observed in 
the low current state does not represent typical behavior, often the current fluctuations are 
larger in the high current state. 

Electron spin-nuclear spin feedback is known to exhibit complex dynamics, including 
hysteresis [HI [IT] , multi-stabilities [2"T] . and fast switching between different stable states 
[5J 122] • This suggests that the observed switching and hysteresis might be due to DNP. 
Earlier experiments attempted to extract the degree of nuclear polarization directly from 
the hysteresis curves [10]. In this interpretation, the maximum nuclear polarization is 
simply given by the size of the hysteresis loop. In our case, Fig. 1(c) would present a 
contradiction to such an interpretation. The field range of hysteresis exceeds 3 T, which is a 
few times larger than the effective field corresponding to full polarization for InAs (~ 1.0-1.5 
T, depending on the effective electronic g-factor). This shows that a straightforward analysis 
of the hysteresis curves does not provide an estimate for the maximum nuclear fields in the 
dots. 

We are able to determine the nuclear polarization directly, by probing the effective Zeeman 
splitting of the electrons. To this end, a continuous wave GHz-frequency electric field is 
applied to one of the gates. The oscillating electric field drives transitions between the 
Zeeman-split spin-orbital eigenstates of electrons. When the a.c. frequency matches the 
Larmor precession frequency / in one of the dots, additional transitions within the (1,1) 
manifold are induced and extra current flows through the double dot. In the presence of 
DNP, the Larmor frequency is given by 2irf = \gfiB{B ext + -B/v)|> where g is the effective 
g-factor in the quantum dot, /xg is the Bohr magneton, is the effective nuclear field, and 
we have set U — 1 for convenience. Therefore, a finite will reflect in a shift of /. 
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FIG. 2: (a) Double dot current as a function of magnetic field (device 2, Vd c = 7.5 mV, and e ~ 2.5 
meV). The gray and white background indicates the high and low current states respectively, (b) 
The frequency of the a.c. voltage on a local gate is swept, while the magnetic field is stepped. At 
each point, we measure 50 ms with a.c. excitation on the gate and then 1 s without. The difference 
in the two measured currents is plotted. This procedure also avoids dragging of the spin resonance. 
The inset shows a zoom-in on the spectral line at the transition from the high to low current state. 
Low a.c. power is used to minimize switches between metastable states induced by the driving. 
The horizontal lines at fixed frequency correspond to photon assisted tunneling enhanced by cavity 
modes in the fridge. This is the same device as studied in [3]. 

The EDSR spectroscopy is performed on a second device, the spectrum is shown in Fig. 
2(b) and the high and low current states in Fig. 2(a). In the low current state we fit the 
observed resonance frequency to 2nf = \gfiB-B ex t\, yielding g = 8.7 ± 0.1. Similar g-factors 
were measured in this device in a regime where hysteresis was not observed [3] . The scenario 
proposed in [TDj , in which DNP compensates the external field fully in the low current state 
to a nearly zero total Zeeman energy, can thus again be ruled out. We also note that the 
extent of the low current state is asymmetric with respect to the zero field axis j23]- This 
effect is further explored in the supplementary material. 

A second important observation is that only a single EDSR resonance is observed in the 
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low current state, while in the high current state we see multiple EDSR resonances. In 
strongly coupled quantum dots we do expect multiple resonances corresponding to various 
transitions within the manifold of (0,2) and (1,1) singlet and triplet states [33]. The ob- 
servation of only a single resonance in the low current state is thus surprising. We at least 
expect to see two EDSR lines due to a g-factor difference between the two dots [3]. 

We also detect a 0.3 GHz shift in the primary resonance frequency at the transition 
boundary [Fig. 2(b), inset]. This suggests that the two current states do have different 
nuclear polarizations, the difference being however rather small (at least in one of the dots). 
How could a small change in polarization alter the current through the double dot by an 
order of magnitude? And why is only a single EDSR resonance observed in the low current 
state? In what follows we propose an explanation based on the gradient in the ^-projection 
of the nuclear fields AB Z N over the two dots [8]. We stress that the physics presented 
here is especially relevant for materials with the strong spin-orbit interaction which couples 
polarized triplets to singlets. 

Let us first consider the basics of transport through a spin-blocked (1,1)— )-(0,2) transition 
in the presence of spin-orbit interaction and in the strong inter dot coupling regime. An 
electron enters the left quantum dot and forms one of the four (1,1) states with an electron 
on the right dot. Out of the four states, only the (1,1) spin singlet is tunnel coupled to 
the (0,2) singlet due to spin selection rules [the (0,2) triplet states are at too high energy 
to play a role]. In the absence of a magnetic field, this results in three blocked states and 
therefore very low current, proportional to the residual escape rate out of the blocked states. 
If a finite magnetic field is applied, the spin-orbital T + (l, 1) and XI (1, 1) are split off by the 
Zeeman energy ±Ez = ±gfiBB ext (g being the average ^-factor of the two dots), see Fig. 
3(a). Spin-orbit interaction then effectively allows for spin-nonconserving tunneling and 
couples T±(l, 1) to S(0, 2), characterized by the energy t so [20]. The remaining blocked state 
T (l, 1) forms the bottleneck for transport and the escape rate out of this state determines 
the current. 

A difference AEz in effective Zeeman splittings in the two dots mixes Tq(1, 1) with S(l, 1), 
thereby unblocking the system at finite magnetic field [5], as indicated in Fig. 3(a). AEz 
is contributed to by a difference in the effective ^-factors of the two dots Ag, as well as by 
a nuclear field gradient along the z-axis. Since T (l, 1) constitutes the bottleneck in the 
transport cycle, a change in AB Z N (and thus in AEz) could indeed have a significant effect 
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FIG. 3: (a) Sketch of the spectrum of the five electronic states. The coupling between the different 
states, as well as the decay rate r out of 5(0, 2) are indicated, (b) Simulated current in the presence 
of an external Zeeman gradient of 1 mT (red trace) and 3 mT (blue trace). We used e = r out = 1.5 
meV, t s = 100 ^eV, t so = 25 /ieV, and g = 9. We averaged over 10 5 random nuclear fields taken 
from a normal distribution with <J ({B^ y,z )' 2 ) = i?^- m ' s ' = 1 mT for both dots. 

on the current. We support this statement by performing transport simulations including 
a Zeeman gradient. Fig. 3(b) shows the calculated double dot current versus B ext , for two 
values of AEz- We used a rate equation model that includes the effects of hyperfine and spin- 
orbit interactions |19j . A gradient of only a few mT is sufficient to increase the current by 
almost an order of magnitude. Note that we have previously reported higher-than-expected 
current levels in the strong coupling regime [TH]. We now propose that this higher current 
was due to the Zeeman gradient over the two dots which was not included in the model at 
the time. 

In InAs double dots, it has been observed that typically Ag/g = 1-10 % [3j [H]. At 
applied fields of 100 mT this mismatch would induce a Zeeman gradient of several mT and 
thus lead to a considerable increase of current due to lifting of the spin blockade [Fig. 3(b)]. 
Low current could arise when a small nuclear field gradient exactly compensates the Zeeman 
gradient due to Ag. The low current state thus has AEz ~ 0. This idea also explains that 
only a single EDSR resonance is observed in the low-current state [Fig. 2(b)] in contrast to 
Ref. [3] where in the same device two resonances corresponding to two dots were resolved. 

In what follows we show theoretically that electron-nuclear spin flip-flops indeed can drive 
the double dot towards AEz — and keep it there. In Fig. 4(a) we plot the spectrum of 
the (1,1) states as a function of AEz close to AEz — 0. The thickness of the lines in the 
spectrum corresponds to the occupation probabilities of the four states one finds when taking 
into account the coupling to the decaying (0, 2) singlet: at AEz = the system has one 
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blocked state \Tq) in which it spends all its time. When AEz deviates from zero, \Tq) and 
\S) acquire a spin-orbital "fb- and .J/b-character and the blockade is lifted. Correspondingly, 
the double dot current increases as AEz moves away from zero. 

At AEz = positive and negative pumping rates are balanced. Hyperfme induced 
electron-nuclear spin flip-flops can cause transitions from the only occupied state \Tq) to 
|tt) an d to with equal probabilities. Since |T ) is an equal superposition of |bl) and 
j these processes will not lead to a net pumping of the nuclear fields. 

As soon as AEz deviates from zero, DNP tries to return the system to AEz — 0. For 
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FIG. 4: (a) Thick blue lines, left scale: Spectrum of the (1,1) states as a function of AEz for 
Ez = —20 /ueV, t s = 60 /ieV, t so = 40 //eV, and e = r out = 1.5 meV. The thickness of the lines 
indicates the occupation probability of the corresponding eigenstate (4 /xeV = 1). Brown dashed 
line, right scale: Current through the double dot as a function of AEz for the same parameters. 
Red dashed arrows: Preferred directions of nuclear spin flips close to AEz = 0. (b) d{AP)/dt as a 
function of AP for three different magnetic fields, negative Ez corresponds to a positive field. We 
used the same parameters as at (a), with AI = 0.7 meV, I = |, N = 10 5 , 1/r = 3 • 10~ 10 /ieV, 
and Ag/g = 0.05. Stable (un)polarized states of the nuclear spin ensembles are indicated. 
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example, when AEz is positive, the most strongly occupied state |T ) acquires a i-t-character 
and transitions starting from |T ) have a preferred direction for spin-flips [21]. This yields 
DNP with net rates R~[ and [red dashed arrows in Fig. 4(a)]: the left dot is negatively 
pumped, the right dot positively. Both these rates indeed drive the system back to AEz = 0. 
Transitions from \S) have the same asymmetry but with opposite sign and counteract this 
pumping. However, since the occupation probability of \S) is much smaller than that of 
\T ), transitions from \T ) dominate. 

This intuitive picture is confirmed by an explicit calculation of all allowed hyperfine 
flip-flop rates. The derivation involves a few straightforward steps, following the approach 
of previous works P, |2T] (see the supplementary material for details). Assuming nuclear 
spin 1/2 for simplicity, we calculate each separate transition rate using Fermi's golden rule. 
Summing over all transition rates, we arrive at an equation of motion for the polarization 
gradient over the dots: 

where AP = \{Pl — Pr), with the nuclear polarization in left and right dots — 1 < Pl(r) < 1- 
N is the number of nuclei in each dot (see supplementary material for a discussion of 
asymmetric dot sizes). A is the average hyperfine coupling energy (AI ~ 0.7 meV for InAs, 
with / the average total nuclear spin). The angle 9 is defined by tan# = (AE z )(Tl ut + 
4e 2 )/2e(t 2 s + t 2 Q ), where AE Z = (Ag)fj, B B ext + AI(AP), and the dimensionless function 



f{9) ~ 1 is given in the supplementary material. In deriving (25), we assumed that \Ez\ S> 
T a for all (1,1) states. We also added a phenomenological nuclear spin relaxation rate 
1/r ~ 0.1-1 Hz. The average polarization P = \{Pl + Pr) is not pumped: DNP merely 
enhances the relaxation rate of P. 



In Fig. 4(b) we plot the pumping curve (25) for AP with realistic parameters, for different 
magnetic fields. We see that at low fields the system has one single stable state close to the 
low current point AEz = (blue curve, Ez = —0.4 meV). At intermediate fields, the high 
current unpolarized state with AP is stable as well (green curve, Ez = —0.7 meV). This 
bistability can manifest itself in switching and hysteresis. At very high fields, DNP becomes 
too weak to counteract nuclear spin relaxation and only the high current unpolarized state is 
stable (red curve, Ez — —1.1 meV). Using parameters from [19] we estimate the maximum 
field for which DNP can stabilize the AP^ « to be 6 T, which is indeed consistent with 
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Fig. 1(c) (see supplementary information for details). 

Often, states of nuclear polarization stabilized by DNP exhibit significantly reduced fluc- 
tuations. Using the spin-flip rates found from the Fermi's golden rule we evaluate the 
variance of the gradient distribution in the polarized state as compared to the unpolarized 
state. We find a relative suppression of the mean-square fluctuations of about ~ 4 • 10 -3 
for the parameters of the blue curve in Fig. 4(d), and of ~ 9 • 10 -3 for those of the green 
curve (see supplementary material for details). In the context of two-electron singlet-triplet 
qubits, this narrowing could lead to an enhancement of the dephasing time T| by more than 
an order of magnitude. 

We thank M. Rudner, H. Bluhm, E. Rashba, Yu. Nazarov, L. Levitov, L. Vandersypen, 
and K. Nowack. This work has been supported by NWO/FOM (Netherlands Organiza- 
tion for Scientific Research), by the Alexander von Humboldt Foundation, and through the 
DARPA program QUEST. 

SUPPLEMENTARY INFORMATION 
THEORETICAL CONSIDERATIONS 

As argued in the main text, the data presented in Figs. 1 and 2 of the main text suggest 
an explanation of the high and low current states in terms of small changes in the gradient 
AB^ of the effective nuclear fields in the two dots. The scenario proposed was that in the 
low current state the nuclear field gradient compensates the difference in Zeeman splittings 
due to the g-factor mismatch of the two dots. If this compensation is exact and the effective 
Zeeman splitting for the electrons in the dots is equal, AEz = 0, then there is one fully 
blocked unpolarized spin-orbital (1,1) state and current is low. It was suggested that a DNP 
feedback mechanism stabilizes this state, explaining the large range of magnetic fields in 
which low current could be observed. A finite effective gradient mixes the two unpolarized 
(1,1) states and thus unblocks the system. A high current state could therefore correspond 
to a situation where the nuclear fields are unpolarized and thus AE Z ~ \(Ag)fj, B B ext \. 

In this Section of the supplementary material, we present the details of our derivation 
of Eq. (1) in the main text and we show how we produced Fig. 4. The derivation involves 
a few simple steps, following the straightforward approach of previous related works [HI [HI 
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FIG. 5: Dragging of the low current state depends on the sweep direction. This Figure shows the 
current in the space of detuning and magnetic field. The detuning cut is through a (1,3)— >(2,2) 
transition in device 1 at VJ} C = —9 mV. Green arrows indicate the sweep direction, the variable 
along the orthogonal axis is stepped after each sweep. The low current state (blue) occupies a 
different area in both detuning and field when (a) the field is swept from positive to negative, (b) 
the detuning is swept from high to low, (c) the field is swept from negative to positive, (d) the 
detuning is swept from low to high. Panels (a) and (c) are part of the same scan in which the field 
is swept back and forth. In these panels the system does not leave the low current state for small 
detuning because the field is not swept high enough to trigger a switch. 
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FIG. 6: The "waiting game". Rules: First, the detuning and dc bias are set. Then the field is 
swept from zero to Bq. As soon as Bq is reached the data are recorded as a function of time, (a) 
and (b): Bq = 300 mT, near the field where the system would switch to the high current state in 
a typical field sweep [device 1, (1,3)— >(2,2) transition]. Here the system typically switches to the 
high current state after a time that varies between 1 to 15 minutes. We never observe a gradual 
current increase between the low and the high current states, the change is always abrupt. If we 
rapidly turn the dc bias off and then back on, the current would return in the high current state at 
these settings. In these plots the high current state is characterized by higher fluctuations. Some of 
the switches are due to charge noise [arrow in panel (a)]. The higher frequency current noise may 
originate from multiple metastable current states, (c) and (d): Examples of a different behavior. 
In panel (c) the system switches back and forth between high and low current states (arrows on 
the right side). These data are taken at the (3,1)— )-(2,2) transition in the vicinity of a high current 
state "island", i.e., a region of high current state surrounded by a low current state, similar to the 
one shown in Fig. S4(b), at a detuning of ~ 5 meV and Bq ~ 0.5 T. In panel (d), the field is swept 
from high to low to reach Bq = 80 mT, such that the initial state is the high current state. Panel 
(e) is an example of low current state persistent for hours. These data are taken in the center of 
the field hysteresis loop [device 1, (1,3) — ^(2,2) transition, Bq = 200 mT]. 
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FIG. 7: (a) High resolution scan of the DNP onset area [device 1, (1,3)— >(2,2) transition, VJj c = 
—7 mV]. In this scan, the magnetic field is swept from negative to positive and detuning is stepped 
after each sweep. At fields of ~ 20 mT and —20 mT, the current exhibits local maxima with a 
weak detuning dependence. We interpret this as being the point where the low current state onsets. 
At fields lower than 20 mT the current is low due to the weakened effect of spin-orbit interaction 
and suppressed hyperfine mixing, see main text. In between the two 20 mT peaks, two smaller 
resonances appear to form a funnel shape as they move closer to zero field at higher detuning, (b) 
is a line cut from (a). 
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FIG. 8: Additional examples of DNP at the (3,1)— >(2,2) transition. The detuning is swept from 
low to high, and the field is stepped after each sweep. Panel (a) demonstrates that the maximum 
extent of the low current state in magnetic field can be a non-monotonic function of detuning. 
Panel (b) shows an even lower current state inside the low current region (dashed line). This lower 
current state forms a funnel shape (appears at smaller fields for higher detuning). The current is 
also suppressed where the funnel shape crosses the baseline near zero detuning. The origin of this 
state is unknown, though a T_ — > S crossing is expected to have a similar shape in detuning and 
field. 
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FIG. 9: Temperature evolution of the low current state [device 1, (1,3)— >-(2,2) transition]. The 
shape of the low current state region is only slightly affected by temperature, up to 2 K. Detuning 
is swept from left to right. 
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FIG. 10: (a) The low current state is shown to extend to different magnetic fields for positive and 
negative fields, (b) When the magnetic field is oriented perpendicularly to the nanowire, the DNP 
pattern is nearly symmetric, (c) For fields directed along the nanowire, DNP is suppressed almost 
completely for one field direction. Asymmetric DNP was confirmed in device 2, where the magnetic 
field was fixed at an angle of 45 deg to the nanowire (Fig. 2). Data in Fig. 1 were obtained with 
the magnetic field perpendicular to the nanowire (out-of-plane) . In that sample, pattern of the low 
current state is highly symmetric. In (a)-(c), the detuning is swept and the field is stepped from 
positive to negative for angles 25, 105 and 185 degrees with respect to nanowire. Measurements 
are from device 3, at the (3,1)— )-(2,2) transition, and V s d = 6 mV. The patterns are unchanged if 
the direction of the stepping of the magnetic field is reversed, (d) Angular dependence of the DNP 
asymmetry, characterized by an asymmetry parameter r/ = (B + + B-)/(B+ — B-), where B + and 
B- are the maximum positive and negative fields to which the low current state extends. The ratio 
7} for the (1,1)— t-(0,2) and (3,1)— >(2,2) transitions (circles, left axis), and the current at B = 200 mT 
(triangles, right axis) for a range of field angles are shown. The inset shows device 3 including 
gates 1-5, the source (S) and drain (D) contacts, and the angle 9 giving the orientation of the 
magnetic field. Recent theory suggests that DNP may exhibit an asymmetry due to spin-polarized 
injection into the double dot [23]. In order to account for the apparent violation of time-reversal 
symmetry, the model assumes that spin-polarized injection is allowed even at zero field due to 
spin-orbit interaction in the nanowire leads. The model predicts maximum asymmetry when the 
external field is aligned with spin-injection field. The simplest interpretation of this spin-injection 
field is to identify it with the spin-orbit field in the nanowire. However, the direction of the spin- 
orbit field extracted from spin blockade anisotropy [panel (d)] in device 3 points perpendicular to 

the nanowire, i.e., orthogonal to the DNP asymmetry. 
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[21]. We treat the hyperfine interaction between the electron and nuclear spins as a small 
perturbation, and we thus first focus on the electron dynamics in the spin blockade setup 
while ignoring the hyperfine interaction with the nuclear spins. We derive a description of 
the effective electron dynamics in the (l,l)-subspace, and from this we extract the quasi- 
stationary occupation probabilities and effective decay rates of the four (1,1) states, as well 
as the current through the double dot. Hyperfine induced transitions between the four 
(1,1) states involving electron-nuclear spin flip-flops are then evaluated with Fermi's golden 
rule. All allowed hyperfine transition rates are combined into an equation of motion for 
the nuclear polarization in the two dots. Finally, we investigate two experimentally relevant 
quantities: (i) the typical fluctuations of the nuclear fields around the polarized stable points 
are estimated from the equilibrium solution of a Fokker-Planck equation for the probability 
distribution of the nuclear field configurations, and (ii) we investigate the predictions of our 
model for the maximally achievable polarization gradient. 

Electron dynamics 

We begin by writing the Hamiltonian for the electrons involved in transport. We start 
with the 5-dimensional basis {\T + ) , \S) , \T Q ) , |T_) , IS02)}, where the first four states are 
spin-orbital (1, 1) charge states and the fifth state 15*02) is the (0,2) spin singlet state. An 
externally applied magnetic field splits the energies of the (1, 1) states. The two spin- 
orbital triplets |T + ) and |TL) will acquire a (renormalized) Zeeman energy of ±Ez- The 
"unpolarized" spin-orbital states \S) and |T ) are coupled by the Zeeman gradient AEz 
over the dots. The total Zeeman term in the Hamiltonian thus reads 

H z = E z { \T + ) (T+\ - |T_) (T_| } + AE Z { \S) (T | + |T ) (S\ }. (2) 

The coupling of the (1, 1) and (0, 2) spin singlets reads in the spin-orbital basis 

H t = -it- \T + ) (S 02 \ + t s \S) (S 02 \ + it z |T ) (S 02 \ + it + |T_) (S 02 \ + H.c, (3) 

where = -^{t x ± it y ), with t x and t y being the coupling coefficients of IS02) to \T X ) 
and \T y ), the unpolarized spin-orbital triplet states along the x- and y-axis respectively 
|2Uj . The "spin-nonconserving" tunneling elements are typically of the same magnitude 
tx, y ,z ~ £ so ~ Oit s , where a can be of the order 1 for materials with strong spin-orbit coupling 
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20] . Finally, we describe the detuning between |S) and |S ( 



02/ 



- £ l'S'02) (S02I 



(4) 



where e is taken positive. In one view, the total Hamiltonian thus looks like 

/ E z -it- \ 
AE Z t s 
H = H z + H t + H £ = AE Z it z 



(5) 



-E z it + 

y it + t s — «t 2 — zt~ — e y 

To this Hamiltonian we apply a basis transformation of the two unpolarized spin-orbital 

states {\S) , \T Q )} -)> {|c) , |6)} as 



l&> 



zt 2 |^)+t s |T ) 



and |c) 



t s \S)+it z \T ) 



(6) 



The new basis state |c) is coupled to the singlet IS02), whereas \b) is not. In the new basis 
{\T + ) ,\c) ,\b) ,\T_) , \S 2)} the Hamiltonian reads 





( E z 











-it~ \ 










AE Z 





+ n 


H = 





AE Z 






















—E z 


it + 




\it + 


v 7 ^ + n 





-it~ 





(7) 



We now assume that the decay rate of IS'02) to the outgoing lead r out is the largest energy 
scale involved. This implies that we can separate time scales in the electronic dynamics, and 
eliminate IS'02) from the Hamiltonian. We write a 4 x 4 Hamiltonian for the (1, 1) states, 
where the effect of the coupling to the decaying IS'02) is twofold: (i) it gives rise to exchange 
terms in the Hamiltonian (H' t )ij = ATi2T2jS / '(r^ ut + Ae 2 ) (note that we have set h — 1 for 
convenience), with T a2 = (a|-£f i |So2), and (ii) leads to decay of all coupled states with the 
decay rates Tj = 4|T i2 | 2 r out / {Y 2 ut + As 2 ), see also [20]. The effective 4x4 Hamiltonian reads 

/ E z \ 
E t AE Z 
AE Z 
\ -E z J 



#(1,1) 



(8) 
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with 

_ 4s(g + g) 

1 out i ^ fc 

where we assumed that Ez ^> 4at 2 e/(r 2 ut + 4e 2 ). Indeed, in the typical regime where 
bistabilities are observed, the energy scale characterizing the strength of spin-orbit effects, 
~ £ 2 /T 0Ut , is comparable or slightly larger than the equilibrium r.m.s. of the nuclear fields 
-B^r m s ' ~ 1 mT. The Zeeman splitting Ez is typically much larger: tens of mT to several T. 
The four basis states of (pi) have the effective decay rates 



r± =T t = 1 -{T X + T y ), T C = T S + T Z , and T b = 0, (10) 



with 



4T d 2 



With #(i i) and the decay rates given as in (10), we can find the stationary occupation 
probabilities of the four states as a function of the parameters in the Hamiltonian. Diag- 
onalization of #(i,i) yields the set of four eigenstates {\T + ) , |1) , |2) , |T_)}, where the two 
unpolarized states read 

|l) = cos^|c) + sin^|6) (12) 

if o . e\ t s + u z .... if e . e\ -t s +u z .... 

= —= cos — h sm - — . t1) H — cos sm - — . It) , 

V2\ 2 2J y/%+%" H V2\ 2 2) y^Tt 2 ^' 7 

|2) = cos °-\b)- sin^|c) (13) 

if e . e\ t s + it z .... if 9 . o\ t s — u z 

= —= cos sm - — . Tl) H t= cos — h sm - — . IT) , 

where is defined by tan# = 2(AE Z )/ E t . The eigenenergies of these states are 

i„ i nz, — tttztt: , „ i„ 1 



E 1 ^-E t + -^Ef + A(AE z y and E 2 = -E t - -\/ E 2 + A(AE Z ) 2 , (14) 
and their effective decay rates read 

2 / n \ 2 



r 1= ^cos-J T c and T 2 = (^sin - J T c . (15) 

We then construct a set of master equations to describe the occupation probabilities p a of 
the four (1,1) states. We assume that the charge transitions (0,2)— >-(0,l)— >-(l,l) take place 
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on a time scale much faster than all decay rates of the (1,1) states, and that refilling of all 
(1,1) is equally likely. We then find the occupation probabilities 

r- 1 

P± = 2rr 1 + [(sin^)-2 + (cos^)-2]r-i' (16) 
Pl 2r t - 1 + [(sin|e)- 2 + (cos|0)- 2 ]r- 1 ' 1 } 

= ^Iggg (18) 

P2 2r t - 1 + [(sin^)- 2 + (cos|0)- 2 ]r- 1 - 1 } 

These occupation probabilities and the eigenenergies found above are the ingredients used 
to plot the spectrum in Fig. 4(a) in the main text (thick blue lines). Under the assumptions 
made above, the current through the double dot simply follows as 

4e 

2r- 1 + [(sini0)- 2 + (cos!0)- 2 ]r-i' 1 ; 

on average 4 electrons are transported in a time which is the sum of the four decay times of 
the (1,1) states. This expression was used to produce the dashed brown plot in Fig. 4(a). 

Dynamic nuclear polarization 

We now include hyperfine interaction into the picture. The nuclear spins in both dots, 
represented by the operators Il,r, are coupled to the electron spins via hyperfine interaction, 

L k L R k R 

(20) 

where the two sums run over all nuclei in the two dots. For simplicity we assume that all 
nuclear spins are equally strongly coupled to the electron spin, which reduces the prefactor 
to the hyperfine coupling energy A divided by an effective number of nuclei Nl,r in the dots. 



The z-components of the coupling in (20) lead to an Overhauser shift of the Zeeman 



energy of the electrons in the two dots. Including this Overhauser shift into our description 
is done by setting E z = gn B B e ^ + \A[{P L ) + (P R )\ and AE Z = (Ag)fi B B cxt + lA[(F L ) - (/£)]. 
In terms of the degree of nuclear polarization of the two dots, — 1 < Pl,r < 1, we can write 
Ez = gfJ,BB ext + AIP and AEz = (Ag)fj,BB ext + AI(AP) introducing the average dot 
polarization P = \{Pl + Pr) and polarization gradient AP = \{Pl — Pr)- 
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The hyperfine flip-flop terms S + I~ and S~I + are the ones that can cause spin exchange 
between the electrons and the nuclei and could lead to DNP. In the main text we explained in 
qualitative terms how the dominant hyperfine flip-flop processes close to the point AEz = 
lead to stabilization of this point. We now have all ingredients at hand to evaluate the DNP 
rates explicitly. We employ Fermi's golden rule to calculate all flip-flop rates, and we add 
for each dot all rates flipping nuclear spin up and all rates flipping nuclear spin down, 

A 2 1 



R d 



TC 



ANd E: 



A 2 

AN d E 2 



riv{r> + |(l| S d |T+)| 2 + I> + |(2| S d \T + )\ 2 

+ r^ 1 |(r_|s d |i)| 2 + r^ 2 |(T_|s d |2)p 
Ptjd {r lP _|<i| sf |r_)| 2 + r 2P _|<2| s+ |r_)| 2 
+ r^ 1 |(T + |5+|i)| 2 + r^ 2 |(T + |5+|2)p 



(21) 



(22) 



where Pf(n,d is the fraction of nuclear spins in dot d which has its spin up(down) [27], and 



we used that E z ^> r 12i< . The occupation probabilities Pi,2,± used in (21) and (22) are the 



ones found in the previous Subsection, (16 )— ( 18 ). These probabilities do not include the 



effect of the hyperfine decay rates themselves, but since we have in the regime of interest at 
most one blocked state, all leading order hyperfine induced corrections the the probabilities 
are of the order Bq/ min{r t , T c }, which we assume to be small. 

We then combine the rates into equations of motion for the polarization in the two dots, 

dP L 2 „ , A 2 



dt 

dP R 
dt 

where we used that 



2 

N~l 
2 

Nr 



( R R ~ R r) 



T 2 C sin 2 9 + AT 2 

IQN 2 E 2 r c sin 2 # + 2ry 
A 2 rg sin 2 + 41^ 
16iV2 E 2 T c sin 2 9 + 2F t ' 



-2P L -sin20), 
-2P K + sin 29), 



(23) 
(24) 



|(1±P^). We finally rearrange (23) and (24) into time-evolution 



equations for the average polarization and polarization gradient, 
d(AP) A 2 r 2 sin 2 # + 4r 2 f Nl + N 2 p 



dt 



32E 2 T c sin 2 9 + 2T t 



NlN 2 R 



N 2 - N 2 
R [2(AP) + sin 29] + ^„M R 2P 



dP 
~dt 



A 2 r 2 sin 2 + 4r 2 f N 2 + N 2 R 



32E 2 r c sin 2 + 2r t 



NlNl 



2P + 



Nr 



NlNl 



NlNl 
N 2 

R [2(AP) + sin 20] 



(25) 
(26) 



In the limit of equal numbers of spinful nuclei in the two dots, Nl = Nr = N, this set of 
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equations reduces to 

d(AP) A 2 r 2 sin 2 fl + 4r 2 /A s n . 1 /A _ 

v - ; = ^ — 5— — 1 2(AP + sin20 - -(AP , 27 

dt 16N 2 E 2 T c sm 2 6 + 2T t [ V ; J r V ; ' V ; 

dP = A 2 rgsin 2 g + 4Tf 1 

dt " 8N 2 Eg T c sin 2 6* + 2r t r ' 1 ' 

where we included a phenomenological relaxation rate 1/r, with r being typically very long 
(often on the scale of seconds). A deviation from Ni = Nr leads to small corrections to the 
pumping curves and stable points. To leading order, one can expect corrections of relative 
magnitude (iV 2 — A^) 2 /(iVf + N R ) 2 , which is typically very small. 



Eq. (2) in the main text is (27), where for simplicity we have set t x = t y = t z = t so . 
Within this approximation, we have for the function f{9) introduced in the main text 

m = jW^M (l + a 2 ) 2 sin 2 g + 4a 4 

+ sin 2 + 2****, (l + a 2 )sin 2 + 2a 2 ' 1 ' 



Eq. (27) is the one we used to produce the plots in Fig. 4(b) of the main text. As already 
can be seen from the plots, for a large range of magnetic fields, the DNP mechanism can 
indeed stabilize the system close to the point with AEz = 0. 

Let us briefly comment here on the differences between the mechanism proposed in this 
work and previously investigated mechanisms. The key difference compared to similar sys- 
tems hosted in GaAs, is the strong spin-orbit coupling, effectively lifting the blockade of the 
Zeeman split-off triplet states. As long as Ez > max{(t 2 /T 0Ut ), AEz} we have at most one 
blocked state [see Fig. 4(a) of the main text]. This allows us to neglect the hyperfine assisted 
escape rates when calculating the pi [the thickness of the lines in Fig. 4(a) of the main text]. 
In a similar situation in GaAs, the hyperfine rates themselves could be the dominant escape 
rates out of the three triplet states, thereby heavily influencing the pi, which could in turn 
significantly affect the DNP [22] • Other groups used elaborate gate pulsing schemes to drive 
one specific hyperfine transition, usually isolating the S —> X+ transition [161 EH]. In this 
case, a small imbalance of the and J/N-components in the electronic ground state can also 
cause a significant pumping of the nuclear field gradient, ultimately leading to a narrowing 
of the distribution function of the nuclear field gradient around AP = [21]. This narrow- 
ing is accompanied by an overall drift of both nuclear polarizations to lower values since 
all nuclear spin flips have the same preferred direction. Apart from this, when driving the 
S — > T + transition, a small asymmetry in the number of nuclei in the two dots (Nl ^ Nr) 



22 



can also have a dramatic effect: the resulting imbalance in the spin-flip rates could lead to 
continuous increasing of the field gradient. In our setup, a small asymmetry between Nl 
and N R only causes small corrections to the pumping equations (25) and (26). 



Maximum gradient achievable 



A useful quantity to extract from Eq. (27) is the maximally achievable polarization gradient 
(AP) max . Using the fact that the gradient at the stable point exactly cancels the field 
gradient caused by the g-factor mismatch, we can relate this maximum as AT(AP) max = 
— (Ag)fiBB max , where P max is the maximum field for which a bistability could be observed. 
We can find this maximum from zooming in on the peak-dip structure. We see that the 
top [28J of the DNP peak corresponds with sin26 l « —1, so that 9 ~ — 7r/4. This makes 
sin 2 # w 1/2, so pumping at this local maximum reads 



d(AP) 



dt 



A 2 rti n 2 t z r l 



4iV 2 P 2 T 2 + 4e 2 rtlr 



{-1 + 2(AP)}--(AP), (30) 



where we introduced the coefficient 



r _ f ( !\1 ^ (3^ 2 + l) 2 t\r 

J V 47 N 2 T ~ 5a 2 + 1 iv 2 r ' 1 ; 

for convenience of notation. When this local maximum in the pumping curve is exactly 
zero, we know that the corresponding polarization gradient (AP) max is the maximum AP 
achievable, and from this we find B max , the largest field for which in principle bistabilities 
could be observed. Setting (AP) max = yields 



\9 



IA 3 



ImnB^l ^ 2A2 + 4r[1 + 4(£/ p out)2] { g^ Bmaxr ( 32 ) 
which is a cubic equation and can easily be solved. 



In Fig. 11 we plot the dependence of P max on the detuning e, using the same parameters 



as were used in Fig. 4(b) in the main text. We see that for large detuning the range of 



magnetic field in which hysteresis can be observed is suppressed. From Eq. (32) we find 
that the large-e behavior of the maximal field is P max oc e~ 2 / 3 . Boundaries of the hysteretic 
regime in the (P,e)-plane with similar shape have indeed been observed in experiment, see 
Figs. S1-S6 but also for instance the data presented in [TO] . 
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FIG. 11: The estimated maximum field for which the pumping curve (27) yields a stable point at 
AEz = 0, as a function of the detuning e. We used the same parameters as in Fig. 4b in the main 
text: AI = 0.7 meV, I = 1/2, N = 10 5 , r out = 1.5 meV, t s = 60 /ieV, t so = 40 /ieV, {Ag)/g = 0.05, 
and 1/r = 3 • 10~ 10 ^xeV. 



Fluctuations of the nuclear fields around the stable points 



To estimate the typical fluctuations of AP around the stable point with AEz = 0, we 



largely follow the approach of [TTJ Since the dynamics found in (27) and (28) for AP 
and P are independent, we can focus on the fluctuations of AP separately. We label all 
possible configurations of the nuclear spins in the two dots resulting in different AP by 
!(iV +i L — N- } l) — 7;(N +j r — N^r), where N± ; l(r) denotes the number of nuclear spins in 



n 



the left (right) dot with spin up or down. This labeling results in 2N different discrete values 
for n, ranging from — N to N. A nuclear spin flip in either of the dots results in a change of 
n by 1. We now construct a Fokker-Planck equation for the probability distribution V(n), 
based on the master equation 

dV{n) 



Ot 



-V(n)[r+(n) + r_(n)] + V(n - l)T + (n - 1) + V(n + l)r_(n + 1). 



(33) 



Here, V(n) gives probability to find the system in the configuration n, and T±(n) is the rate 
at which the nuclear spin baths flip from the configuration n to n ± 1. Using the fact that 
the number of configurations 2N is large, we can go over to the continuous limit, yielding 



dV d { Id 



(34) 
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which indeed is a Fokker-Planck equation. 

Since the functions r±(n) are smooth on the scale of n, we can use that |d n r±(ra)| <C r^. 
In terms of the polarization gradient AP = n/N this allows us to find the approximate 
solution valid close to the stable point AP = AP , 

V{AP) « V(AP ) exp | - (AP 2 ~f Po) } , (35) 

where the coefficient a 2 is defined as 



a 2 




d r . - r_ 



d(AP) r + + r_ 



-i 



1 r, + r_ 



2N mh*(F--r- 



(36) 

AP 



AP / ~ J d(AP) 

In the last simplification, we used the fact that (r + — r_)| APo = 0. We see that this solution 
for V(AP) is a Gaussian distribution with a width a. 

The rates of change of n can be expressed in terms of the spin flip rates as T + = + R^ 
and T_ = R£ + P^, from which we can write an explicit expression for a. To arrive at a 
rough estimate, we evaluate the fluctuations at the point where exactly 9 = 0. This yields 



a 



2 



J E t [N 2 E 2 (T 2 out + As 2 ) + a 2 A 2 T out t 2 r} 

2N E t [N^E 2 (T 2 out + As 2 ) + aWT^r] + 2a 2 AHY mt t 2 s r 



(37) 



In reality, 9 is not exactly zero at the stable point. Expanding (36) around the point where 
9 = 0, we find to leading order corrections of ~ 5(AP)[(Ag)fj,BB ext /E t ] to both (r + + T_) 
and <9(AP)(r_ — T + ), where 5(AP) is the deviation of AP from the point with 9 = 0. Due 
to the shape of the pumping curve, 9 will always lie within — 7r/4 < 9 < 7r/4, which sets the 
scale for the maximal <5(AP). From this estimate, it follows that both relative corrections 
are of the order < \(Ag)fj,BB ext /AI\, which is indeed small. 
In the regime where 

AI » E > and fAf^ft»? < 38 > 

the fluctuations can be estimated as 

„ a ^ e(t 2 + t 2 J , N 2 E\ TL + 4e 2 



NAI(T 2 out + Ae 2 ) \ + tA* T out tl 



For all experimentally relevant parameters we expect (38) to hold. 



If we combine these estimates for the fluctuations of the nuclear field gradient with the 



simple expression we found for the double dot current (19), we can investigate what our 
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model predicts for the current fluctuations induced by the nuclear field fluctuations. At first 
sight, one might expect lower current noise in a state with small fluctuations. However, 
since the current is a very sensitive function of AEz close to the polarized stable point, 
small fluctuations around 9 = can have a dramatic effect on the current. Expanding the 
current around 9 = up to second order in 5(AP), we can estimate the magnitude of the 
current fluctuations around the polarized stable point as 

A 2 T 2 V 

5I ~ 4e 1 ° U V , (40) 
E t e 



which under the condition of (38) simplifies to 



AIT out / , N 2 E 2 z r 2 out + 4e 2 



" - nvf i 1 + witsrJ ■ (41) 

To put this in perspective, we evaluate the current for AP — > ±oo, i.e., for 9 — > ±7r/2, 

2T ont E t 1 + a 2 



1 max 



e ° ut 1 ' „ . (42) 
e l + 3a 2 V ; 



This allows us to express the relative fluctuations in the current: the magnitude of the low 
current fluctuations scaled to the difference in high and low current, / max - We thus find 

51 AI 1 + 3a + A^i £L + 4^\ _ (4g) 



J max 2iVE 4 1 + a 2 V r A 2 T out t 



'so 



As written above, the typical regime where strong bistabilities are observed is that of 
strong, but not too strong, coupling, where the energy scale characterizing the strength of 
spin-orbit effects in the (l,l)-subspace is comparable or slightly larger than the equilibrium 
r.m.s. value of the nuclear fields. More quantitatively, we expect to be in the regime where 
^so/rout ^ Bj^ m - S - ~ AI/yN. We see that when a ~ 1 (strong spin-orbit coupling), we have 



E t ~ t 2 Q /r out and thus find for the prefactor in (43) 



AI 1 + 3 « 2 <^ (44 ) 



2NE t 1 + a 2 ~ v/jV' 

which is small, typically ~ 10 2 10^ 3 . The actual ratio SI/I max then depends on the second 
term within the brackets, which is ~ 1-10 for the parameters used in the main text. In this 
regime the relative current fluctuations are thus indeed small. 

This ratio however, does not say anything about the relative magnitude of the fluctuations 
around AP = and those around the polarized state. In fact, for the parameters we 
used in the main text, we find that the model predicts the current fluctuations in the low 
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current state to be larger than in the high current state, which is consistent with Fig. 1(c) 
of the main text. Other data sets however (such as shown in Fig. S2) seem to indicate 
opposite behavior. A full understanding of these different observations requires a much 
more quantitative examination of all parameters which goes beyond the scope of this work. 
We just note that the simple model presented here in fact does not exclude the possibility 
of having larger current fluctuations in the high current state. 
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